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In order to understand the gluon contribution to the nucleon spin, 
some experiments can study the production of hadrons at high transverse 
momemtum from lepton-nucleon or nucleon-nucleon scattering. RHIC has 
recently measured such double spin asymmetries All{pt) for pion produc¬ 
tion at high center of mass energies HI 0. and inclusion of its data to 
global hts based on NLO collinear pQCD calculations gives some con¬ 
straints on the gluon polarization in the range 0.05 < Xg < 0.2. [3] 

To complement these results at COMPASS range, we will present 
preliminary COMPASS results on double longitudinal spin asymmetries 
All{pt) for single hadron production measured on deuteron and pro¬ 
ton target at < 1 GeV^, pt > 1 GeV/c and center of mass energy 
y/s ~ 18 GeV. All GOMPASS data taken from 2002 to 2011 by scat¬ 
tering 160 GeV polarized muons on longitudinally polarized ®LiD and 
NH 3 targets have been used, and the number of hadrons collected with 
Pt > I GeV/c for this analysis amounts to about 10 millions. The ob¬ 
tained asymmetries will be compared to theoretical predictions of at NLO 
without gluon resummation calculation. 
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Introduction 


The data used in this analysis come from the COMPASS experiment, located at 
CERN. It uses a polarized muon beam coming from the SPS at 160 GeV (2002- 
2007) or 200 GeV (2011) scattered off a polarized ®LiD (2002-2006) or NH 3 (2007- 
2011) target. The GOMPASS spectrometer allows us to select precisely the data at 
Pt > 1 GeV/c and < 1 GeV/c^ with also other kinematical and geometrical cuts 
to place oneself in the theoretical settings and have cleaner data. Among these cuts, 
the main ones are on y the fractional energy loss (0.1 < y < 0.9) to avoid large 
radiative contributions and region with imperfect Monte Garlo; on 2 ; the fractional 
hadron energy ( 0.2 < 2 ; < 0 . 8 ) to avoid current fragmentation and exclusive events; 
and hnally on 9^ ( 0.01 < 6h < 0.12 rad) which is the only region well covered by the 
spectrometer (especially for 2002-2004 data). 

The previous studies of these data were using the Monte Garlo generator PYTHIA 
to unfold the partonic information at LO from All. This new study allows to ex¬ 
tract these information at NLO thanks to a new theoretical framework that will be 
followingly presented. 


1 Theoretical Framework 


This analysis is based on collinear perturbative QGD calculations done by B. Jager, 
M. Stratmann and W. Vogelsang in 2005 [1]. Unpolarized and polarized cross-sections 
can be computed in function of px and p (the pseudo-rapidity) by convoluting nucleon 
PDFs, with perturbative partonic cross-section, fragmentation function distribution, 
and a last part A/^ corresponding to a 'muon distribution function'. This last dis¬ 
tribution leads to two different parts of the calculations: 


Direct p—contribution illustrated in Fig. 
convolution of Weisacker-Williams 


^with A/^^(a;a,/i/) = ^^P^^,{y)5{l- 

)robability with a Dirac distribution. 


• Resolved p—contribution illustrated in Fig. gwith A/„^(a:„,/i/) = fAP^f,{y)Af2{x^ 
p/): convolution of Weisacker-Williams probability with a perturbative photon- 
parton distribution. 


dAcr^ _ Sa,fe,e A/^ ® Af^ 0 dAag^b^cX Dc _ dAadir + dA^res 

dah ^ J2a,b,c fa ® fb ® daa,b^c,X G) da dir + dares 

The collinear assumption implies that the measured px should arise from hard 
subprocesses which can then only be higher order ones. The leading order a'^as 
includes 2 direct processes: the photon-gluon fusion {^*9 99) which probes the 
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gluon distribution and the QCD Compton ( 7 *g —)■ qg), but also other resolved pro¬ 
cesses taking into account the gluon distribution. The next-to leading order 
with a multitude of direct and resolved processes, makes appear leading logarithms 
aaglog‘^^~^{w/{l — w)) which can become large at low center of mass energy and need 
to be checked out. 

For this purpose a preliminary study has been done on the unpolarized cross- 
section |5] to check if this theory used in RHIC (a/s ~ 200 GeV) was still applicable 
in COMPASS kinematics (a/s ~ ISGeV). The study concludes that in COMPASS 
kinematics, gluon resummation is needed to explain the cross-section and its y — 
dependency (Fig. i 0- 
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Figure 3: Unpolarized cross-section measurements compared with theoretical calcu¬ 
lations at LO, NLO and NLO-I-NLL (left), and ratios experiment/theory for NLO 
and NLO-I-NLL (right) 
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2 Measurement Method 


The asymmetry extracted from COMPASS data is a lepton-nucleon semi-inclusive 
asymmetry based on hadron counts. This asymmetry is diluted by [3 composed of 
3 factors: / the dilution factor of the target, the beam polarization and Px the 
target polarization (divided in cells with opposite polarizations). These different cells 
give us two hadron counts and N^) and two other after rotation of the solenoid 
field {Nui and Nd’)- 


N — Nrl 

Nx = 0aj;nj.(To(l + {f ■ Pfj^ ■ Px)All) ^raw = ^ 

I\u + I\d 

To remove more efficiently the acceptance effects between the cells, a second order 
method is used to extract the asymmetry from the quantity 6: 

. _ Nu- Nd' ^ (1 + {I3u)All){ 1 + {fid') A ll) 

Nd-Nu' {1 + {l3d)ALL){^ + {/3 u')^ll) 

To statistically optimize {An) one weights each hadron by tc = / ■ to get 

^ /(yP- 
-A,* ~ V ■ 


^ _ J2wu ■ J2wd' ^ (1 + {f3u)w ■ All){1 + {f3d')w ■ All) 

J^Wd Wu' (1 + {I3d)w ■ All) (1 + {/3u')w ■ All) 

Asymmetry is hnally extracted for each group of data separated by two solenoid 
field rotations and then averaged to remove long term instability: 


All = 


Si ^Ai Ai 


(4) 


3 Systematic Studies 

A various of qualitative studies have been done on ’/a/se’ asymmetries (non physical 
asymmetries supposed to be 0): 

• Misconhguration asymmetries (grouping of data with identical spin state sup¬ 
posed to give a 0 asymmetry) showing the apparatus asymmetry. 

• Difference of asymmetry between left part and right part of the spectrometer, 
and difference of asymmetry between top part and bottom part of the spec¬ 
trometer, showing efficiency anisotropy of the spectrometer. 
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Difference between upstream cell and downstream cell showing acceptances ef¬ 
fects and polarization inhomogeneity of the target. 


• Difference between day and night asymmetries showing effects of thermal ex¬ 
pansion and electronic noise fluctuation. 

In conclusion, no false asymmetry has been detected for this analysis under the sta¬ 
tistical uncertainties. 

Other quantitative studies have been pursued: 

• A multiplicative error study: + (y) ~ Q.Ql All 

(where Pt = {Px)x)- 

• A time stability study through pulls distribution Ar = , from which one 

can extract a systematic error value of the same order as the statistical one for 
this analysis. 


4 Results 

Asymmetries have been calculated for each year to hnd a good compatibility between 
different setups for deuteron and proton target separately (Fig. |^. 




Figure 4: for 2002-2006 COMPASS deuteron data (left) and for 2007-2011 

COMPASS proton data (right) 

These results are on the whole compatible with 0 except for high px deuteron and 
low Pt proton. This last observation was in some way expected and observed in other 
analyses. 
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5 Comparison with Theoretical Calculations 

These results have been compared with theoretical calculations using the same code 
as the one in 2005 with some exceptions: 

• Unpolarized PDFs come from CTEQ6 (as previous calculations). 

• Fragmentation functions distribusions come from DSS. 

• Unpolarized photon-parton distributions come from GRS and polarized ones 
from W. Vogelsang code. 

• 3 different sets of polarized PDFs have been used to show the effect of AG on 
All{pt)'- 2 extreme sets of GRSV [7] and DSSV++ 2014 |3]. 

These calculations have been of course done to match kinematics of GOMPASS, 
i.e. a Weisacker-Williams probability with Q'^^x — ^ {GeV/cY and cross-sections cut 
with 0.1 < 2 / < 0.9 and 0.2 < 2 ; < 0.8, and integrated over rjcms G [—0.1, 2.4] which 
corresponds to the 6h cut of the data. 




Figure 5: Gomparison between experimental An and NLO theoretical asymmetries 
for different sets of polarized PDFs on the left and right) 

These hrst comparisons (Fig. are only preliminary, knowing that the gluon 
resummation should change the theoretical calculations. However the deuteron data 
points would indicate a large positive AG, whereas the proton data points cannot be 
explained by any set of PDFs, how large the AG is. This disagreement between the 
different target material results makes difficult to interpret them. 


Prospects 

Even though this data analysis is already well advanced, some work remains to do to 
get clearer data by separating charges and flavors. For this purpose it is foreseen to 
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extract from the same data Al'l to compare with the theoretical 

calculations. Getting and A’ll is also possible but will bring us less information 
due to the higher statistical uncertainty of these data. 

Of course this comparison with theoretical calculations is just a hrst step to il¬ 
lustrate the link between An and AG, but these results should be in the future 
integrated to world data hts in order to reduce uncertainties over Ag{xg) in COM¬ 
PASS Xg region which is not yet well determined for these data. 

Even though this data analysis is already well advanced, some work remains to do 
to get clearer data by separating charges and flavors. For this purpose it is foreseen to 
extract from the same data A^l, An to compare with the theoretical 

calculations. Getting A|jj and A|jj is also possible but will bring us less information 
due to the higher statistical uncertainty of these data. 

Of course this comparison with theoretical calculations is just a hrst step to il¬ 
lustrate the link between An and AG, but these results should be in the future 
integrated to world data hts in order to reduce uncertainties over Ag{xg) in COM¬ 
PASS Xg region which is not yet well determined for these data. 
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